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Abstract 

We study the potential of a static quark anti-quark pair in the confinement "phase" 
of the SU(2) Higgs model. Around separation r&, the confining string of the gauge field 
breaks by formation of a dynamical pair of light quarks. The ground state and first 
excited state static potentials are determined by a variational technique from a matrix 
correlation in which suitably smeared gauge and Higgs fields enter. Our results at 
(3 = 2.4 clearly show string breaking (r& ~ 1.9ro). The investigation of properly defined 
overlaps confirms the interpretation of string breaking as a level crossing phenomenon 
between string-type and meson-type states. We study the scaling properties of the static 
potentials along a line of constant physics, varying the lattice spacing by a factor of 2. 
Our results show compatibility with scaling within tiny errors. 



COLO-HEP-444 
DESY 00-078 
May 2000 



' e-mail: knechtli@pizero.colorado.edu 
^ e-mail: sommer@ifh.dc 



1 Introduction 



The potential of a static quarkP] anti-quark pair in a Yang-Mills gauge theory has been 
computed by Monte Carlo simulations on the lattice. The observable from which the 
static potential is extracted is the Wilson loop. The results for gauge group SU(2) 
and SU(3) HQ, close to the continuum limit, show a linearly rising confinement 
potential at large separations of the static charges. When the Yang-Mills gauge theories 
are coupled to matter fields in the fundamental representation of the gauge group, the 
static potential is expected to flatten at large distances: the ground state of the system 
is better interpreted in terms of two weakly interacting static-light mesons which are 
bound states of a static and a dynamical quark. The dynamical quarks are pair-created 
in the strong gauge field binding the static quarks. This phenomenon is called string 
breaking or screening of the static charges. The name "string" refers to the gauge field 
configuration which confines the static quarks and leads to the linear confinement in 
pure gauge theories. 

In recent attempts in QCD with two flavors of dynamical quarks PJ^,@J§,g. 1C , 11 ] , 



the flattening in the static potential determined from the Wilson loops was not visible. 
The string breaking distance r^ around which the static potential should start flattening 



off, could nevertheless be estimated in the quenched approximation of QCD to be 12 . 13 ] 



r b « 2.7r , (1) 

where ro ~ 0.5 fm is a scale conveniently defined from the force, F(r), between static 
quarks Q: 

F(r )r 2 Q = 1.65r , (2) 

The QCD results so far show a linear rise of the static potential for distances beyond rb- 
It has been commonly argued that the problem is the poor overlap of the Wilson loops 
with the ground state of the system. The investigation of the static potential in models 
other than QCD is therefore relevant in order to understand its origin and identify the 
reason for the failure of the method used to extract it in full QCD. 

First studies of string breaking were performed with a hopping-parameter expansion 
in SU(2) gauge theory with Wilson fermions [jlj]. In the Schwinger model (QED 2 ), the 
exact solution for the static potential can be given in the limit of zero fermion mass 
fl5|| : V(r) = (ey / 7r/2){l — exp( — er/y 7 ^)}, where e is the charge of the static sources. 
String breaking was established by numerical simulation in the Schwinger model p^Jl7[| . 
Numerical evidence of the screening of the static potential was also found in the U(l) 
Higgs model (scalar QED) in two dimensions [^] . The flattening of the static potential 
at large distances is also expected in the confinement "phase" of the SU(2) Higgs model. 
Indeed, early simulations yielded some qualitative evidence for string breaking |19|,2C|. 



A static quark (anti-quark) is an external source in the (complex conjugate of the) fundamental 
representation of the gauge group. 
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String breaking can also be studied in Yang- Mills theories using static sources in the 
adjoint representation of the gauge group. The gauge field itself is responsible for the 
screening of the sources and the formation of hadrons called "gluelumps" . An important 
numerical investigation concerning this screening has been carried out by C. Michael 



in 1 21], where it has been noted that string breaking can be a mixing phenomenon. 
The static potential is extracted from a matrix correlation in which two types of states 
enter, the adjoint string and the "two-gluelump" . However, due to large errors, no 
clear evidence for string breaking could be given. The first numerical evidence for 
string breaking in non-Abelian gauge theories with dynamical matter fields was given 
using the mixing method in the four-dimensional (2^] and three-dimensional [^3| SU(2) 
Higgs model by the computation of the potential between static quarks. Most recently, 
the extraction of the static adjoint potential in the three-dimensional [^,^] and four- 
dimensional [pf| SU(2) Yang-Mills theory shows also evidence for string breaking. 



In full QCD, string breaking has been seen at finite temperature [27], where the 
static potential can be extracted from Polyakov loop correlators. A recent investigation 
in zero-temperature QCD with two flavors of dynamical quarks [28] gives some indica- 



tion that string breaking can be observed in the potential determined from a matrix 
correlation containing string-type and meson-type states. The main problems are the 
computational costs of the light quark propagators entering the matrix correlation: the 
maximal variance reduction method of [2S] has been applied in [28]. 



In this article, we present results for the spectrum of static-light mesons and for 
the static potential in the confinement "phase" of the SU(2) Higgs model. The method 



of computation of the static potential is the same as in our first work [22] but it is 
explained here in much more detail. In Sect. g we present a detailed investigation 
of the spectrum of the static- light mesons. This is relevant not only for determining 
the asymptotic value of the static potential, but also for finding a suitable smearing 
procedure for the Higgs field to be used in the computation of the static potential. In 
Sect. ||, we show our results for the ground state and the first excited static potential 
obtained at (3 = 2.4, with a spatial resolution two times better than at (3 = 2.2 used in 
|2l 



In order to study the interpretation of string breaking as a mixing phenomenon [3C] 
we first properly define overlaps of the string and two-meson states with the ground and 
first excited energy eigenstates. Their dependence on the distance r then establishes 
string breaking as a level crossing phenomenon. In Sect. [3], we study the scaling of the 
static potentials by comparing the results at (3 = 2.4 with (3 = 2.2 on a line of constant 



physics ]31|]. The dependence of the static potentials on the value of the Higgs quartic 



self-coupling is also investigated. 

2 The spectrum of the static-light mesons 

The investigation of the spectrum of static-light mesons is an important study to be done 
before the extraction of the static potential. The system composed by a static quark 
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anti-quark pair is expected to be described at large separation of the static sources by 
two weakly interacting static-light mesons, which are bound states of a static quark 
and the dynamical matter field. Denoting by [i the mass of the lowest meson state, the 
static potential Vo(r) is expected to approach the value 

lim V (r) = 2//. (3) 

r — >oo 

Therefore, the mass /U basically determines the string breaking distance around which 
the potential starts flattening out. 

The extraction of the static-light meson spectrum is representative for the vari- 
ational method that we employ also for the extraction of the static potential. We 
constructed a large basis of operators that create one-meson "states" when applied to 
the vacuum. The variational approach chooses the best linear combinations of these 
"states" which approximate the energy eigenstates. Through the determination of the 
meson spectrum we can therefore gain information about a suitable way of constructing 
the meson-type states. Then, we use this information in Sect. |3] to construct a basis of 
states for the determination of the static potentials. 



2.1 The matrix correlation 



In [32 1, an Hamiltonian formalism for the SU(2) Higgs model is constructed. Along the 
lines of [^] , a transfer matrix operator is defined and its strict positivity for k > and 
A > is provedP]. This property is equivalent to the reality of the energy spectrum in any 
sector of the Hilbert space. Different charge sectors of the Hilbert space are defined by 
the transformation property of the states under gauge transformation. Through Gauss' 
law, this gauge transformation property is related to the presence of static charges in 
some irreducible representation of the gauge group. 

The static-light meson states belong to the charge sector with one static charge in 
the fundamental representation of the gauge group localised at a certain space position 
x. The meson states (i = 1,2,3,...), are described by (composite) fields Of l (x), 
which are constructed with field variables taken at equal time xq and transform under 
a gauge transformation {A(x) G SU(2)} according to 

[Of' A (x)] a = Al a ,(x)[Of(x)] al , (4) 

where a, a' = 1, 2 are color indices. The obvious choice is to take O m (x) to be the Higgs 
field $(x). We are also going to consider non-local linear combinations which take into 
account contributions from Higgs fields at neighboring sites (smeared fields) and more 

2 One can show that the partition function of the SU(2) Higgs model satisfies the property 
Z(/3,k,\) = Z(P, — K, A). There is a mapping of the observables such that the expectation values 
at positive ft are reproduced by expectation values at negative ft. This motivates the restriction of the 
parameter region to the values k > 0. For A = 0, strict positivity of the transfer matrix holds for 
< ft < 1/6. 
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general composite fields, with the intent to model the true wave function of the meson. 
From the basis of meson-type fields Of l (x) a matrix correlation 

cff(t) = ([oy(x + t6)XU(x,x + t6)l[o™(x)} b ). (5) 

is constructed representing the transition amplitude over a time interval t from the 
meson state i to the meson state j. The static charge is represented by a straight 
time-like Wilson line U(x,x + tO)^ connecting x with x + tO. With the help of the 
reconstruction theorem proved in [p2| , it is possible to show that in the limit of an 
infinite physical time extension T of the lattice^] the correlation matrix eq. (||) can be 
written like 

Cff(t) = ^2(j\a)(a\i)e- tw <*, (6) 

a 

where \a), (a = 0, 1, 2, ...) are the orthonormal meson energy-eigenstates with energies^ 
W a , W a < W a+ \. The matrix correlation eq. (||) can be measured in a Monte Carlo 
simulation: now, we describe the variational method for extracting the meson energy 
spectrum from it. 



2.2 Variational method 

For matrices of the type in eq. (^) a general lemma for the extraction of the energies 
W a has been proved in [34]. In this reference, a variational method is proposed, which 
is superior to a straightforward application of the lemma. It consists in solving the 
generalised eigenvalue problem: 

y] Cjj (t)v a j (t, tp) = X a (t, to) ^ Cij(to)v a j(t, to) , X a > X a+ i , (7) 

j j 

where to is fixed and small (in practice we use to = 0). The generalised eigenvalues 
X a (t, to) are computed as the eigenvalues of C = C(to)~ l l 2 C{t)C(to)~ l l 2 and the vectors 

Va,i = yy3{to) l/2 )ijV a j{t,to) With ^V^iVa'j = 5 aa i (8) 
j i 

are the orthonormal eigenvectors of C. The positivity of the transfer matrix ensures 



that C(t) is positive definite for all t. In [34] it is proven that the energies W a are given 
by the expressions 

aW a = \n(X a (t-a : t )/X a (t,t )) + o(e- tAW «) , (9) 

where AW a = min|W a — Wn\. It is expected that, for a good basis of states, the 
coefficients of the higher exponential corrections in eq. (^) are suppressed so that the 

3 In practice, the limit T — > oo is reached when Tmu 3> 1, where mu is the Higgs mass defined as 
the mass gap in the zero charge (gauge invariant) sector of the Hilbert space. 

4 We normalise the vacuum energy to be 0. 
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energies can be read off at moderateiy iarge values of t from the right-hand side of 
eq. (|). 

The variational method eq. (|7]) and eq. @ is our standard method for extracting 
the energy spectrum. What we have stated here about this method is valid for any 
charge sector of the Hilbert space. One has to start from a basis \i) of states belonging 
to that charge sector. The matrix correlation Cij(t) corresponds to matrix elements 
(j\TF n \i), n = t/a, of powers of the transfer matrix operator TT appropriate for the 
charge sector. In [p^] , this correspondence is derived in detail for the sector with a static 
quark anti-quark pair: the energy spectrum are the static potential and its excitations. 



2.3 Smeared Higgs fields 

We studied different bases of meson- type fields Of [ (x) by measuring in Monte Carlo 
simulations the matrix correlation function Cij(t) defined in eq. (|B|) and computing 
from it the energy spectrum of the static-light mesons using the variational method 



described in Sect. 2.2. Our aim was to find the best field basis for describing the ground 
state of the static-light mesons. For these studies we simulated the SU(2) Higgs model 
on a 20 4 lattice with parameters j3 = 2.2, k = 0.274 and A = 0.5. This parameter point 
is in the confinement "phase" of the model and is the point that we used in our first 



work [22]. The measurement of the matrix correlation is improved by the use of the 



one-link integral method [35] 



We first studied a basis containing the fundamental Higgs field <fr(x) and smeared 
Higgs fields obtained by iterating the application of a smearing operator Si to the Higgs 
field. The smearing operator S\ is defined as 

Si*(x) = $(x)+ ]T U(x,y)$(y), (10) 



where U(x,y) is the link connecting y with x. Iterating the smearing operator Si we 
obtain smeared Higgs fields $^ m ^(x) = S^^ix), (m = 0, 1,2, ...), with different smear- 
ing levels m (m = corresponds to the fundamental Higgs field in the Lagrangian). 
We measured a matrix correlation function with a basis of smeared Higgs fields corre- 
sponding to smearing levels 0,1 and 2 of Si. The result for the ground state extracted 
according to eq. (||) is shown in Fig. [I] (triangles). We were not able to reach a plateau 
for the ratio ln(X a (t — a)/X a (t)) within the range of t considered (up to 8 in lattice 
unit). 

Then, we investigated a larger basis of meson-type fields, defining in particular a 
smearing operator S2 as 

5 2 $(x) = V{V<5>{x)+V U(x,y)$(y) + 

\x— y\=V2a 

x o=yo 

V U(x,y)<f>(y)} , (11) 

\x — y\ — y/3a 
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Figure l: Here, we compare the extraction of the mass ji of a static-light meson using 



different smearing operators defined in eq. (10) and eq. (11). The smearing levels used 
are indicated by n and m. 



where = 3>/V and U(x,y) represents the average over the shortest link con- 
nections between y and x. Through iteration of 52 we obtain the smeared Higgs fields 
$2"^ 0*0 = S 1 ™ 3>(x), (m = 0, 1, 2, ...). We considered the following basis of meson-type 
fields Of{x), i = 1,2, ...,11: 



Of{x) 
Of{x) 

Of{x) 

Of{x) 

Of{x) 
Of{x) 



V$(x) , 

V £ U(x,yMy), 



\x — y\—a 
X0 = V0 



V ]T U(x,y)*(y), 

| x — y | = V2 a 

x o=vo 

V U(x,y)*(y), 

\b— y\=V3a 
^0=V0 

^~ 4) (x), 1 = 5,6,7,8, 
1 3 

$(x) x - (x - ak)U(x - ak, k)$>(x) + 



(12) 
(13) 

(14) 
(15) 
(16) 
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Figure 2: Here, we compare the extraction of the mass [i of a static-light meson using all 
the fields eq. (|l~2|)-eq. (19) and only the smeared Higgs fields corresponding to smearing 
levels 1,2,3,4 of the smearing operator 52 defined in eq. dill). 



® J (x)U(x, k)$(x + ak)} 



or (x) 



Of x {x) 



^)x- J2 {Pki{x) + P k i{x - ak) + 

l<fc<K3 

P k i(x -ak- at) + P k i(x - at)} , 
$(x) x ($ t (x)<I>(x)) , 



(17) 



(18) 
(19) 



where in eq. (|l|) Pja(x) = U(x,k)U(x + ak,l)U\x + al,k)U^(x,l) and we use the same 
notation conventions as in p2[| . In Fig. ||[ the result for the extraction of the mass of 
a static- light meson using the fields 0^(x), i = 1, 11 is shown (triangles). Note the 
enlarged scale on the y-axis as compared to Fig. |]. We obtain a nice plateau already 
at moderately large values of t. The situation remains practically unchanged (also 
the statistical errors) if we remove from the basis all fields except the smeared fields 
obtained by iterations of the smearing operator S^- This means that this smearing 
procedure contains all relevant features for describing the ground state which could be 
obtained by using the larger basis. 

When the generalised eigenvalue problem eq. (Q) is solved, the optimal linear com- 
bination of the basis fields Of^(x) describing the ground state can be expressed in terms 
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Figure 3: Here, we show the extraction of the mass of the ground and first excited 
meson state at j3 = 2.4. The basis of meson-type fields was obtained using the smearing 
procedure S2 eq. (|ll|) with smearing levels m = 1,3,5,7,10,15. The simulation was 
performed on a 32 4 lattice and the statistics is of 800 measurements. 

of the components of the vector vq as v o,iO^(x). Therefore, we call vq the ground 
state wave function. An interesting fact we can learn from vq, is that the field Oq}, with 
nearest neighbor contributions, has a very small coefficient vq^- This explains our orig- 
inal difficulties in extracting the meson ground state. In Fig. [l], a direct comparison of 
the smearing operators S\ and S2 shows clearly that the contributions from the excited 
states are much more suppressed when we use S2 (circles) . 



2.4 The meson spectrum at (3 = 2.4 

In Fig. |, we show the results for the static-light meson spectrum that we obtained for 
the parameters (3 = 2.4, k = 0.2759, A = 0.7 (in the confinement "phase") on a 32 4 
lattice. More details about this simulation will be given in Sect. ||. For the measurement 
of the matrix correlation function we used a basis with the six fields 

$^ m) (x) , m = 1, 3, 5, 7, 10, 15 , (20) 

obtained by iterating the smearing procedure S2 in eq. (|Tl|) . As we will see in Sect. ||, 
the lattice spacing at (3 = 2.4 is reduced by almost a factor two with respect to the 
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Figure 4: Here, we show the same data as in Fig. |3| but plotted against the correction 
term exp(— tAW) in eq. @. We use AW = [i* — [i for the ground state and AW = 
fj,** — ji* for the first excited state, see (21). 



lattice spacing at /3 = 2.2. Therefore, at /3 = 2.4 smeared fields with high smearing 
levels m are expected to play a more important role than at f5 = 2.2. This expectation 
is confirmed by the simulation. In order to determine with confidence the static-light 
meson masses, we plot in Fig. Qthe logarithmic ratios on the right-hand side of eq. (||) as 
functions of the correction terms exp(— tAW). This enables us to choose the best time 
t for reading off the masses from the logarithmic ratios and to estimate the systematic 
errors associated with this choice. For the mass of the ground state, we must take the 
largest value t/a = 9. For the mass of the first excited state, we can take t/a = 8. In 
both cases, the systematic errors^ are of the same magnitude as the statistical errors. 
However, these errors are small. The results for the meson spectrum are[] 

a^i = 0.517(2), a/x* = 0.88(3), afi** = 1.21(9). (21) 

We note that the convergence of the right-hand side of eq. @ is not so "critical" in the 
case of the static potentials considered in Sect. ||. 

5 The systematic errors for the masses are estimated from the difference between the mass read off 
at the chosen value of t and the crossing point of the dotted lines in Fig. ^ with the y-axis (t = oo). 

6 We use the notation fj,, fj," and fi** for Wo, W\ and W2 respectively. 
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3 String breaking and mixing 



As mentioned in the introduction, the basic point concerning the determination of the 
static potential has been first noted by C. Michael [21] in a study of the SU(2) static 



adjoint potential. The energy-eigenstates of the system composed of a static quark anti- 
quark pair and of light dynamical matter fields are well described by a superposition 
of string-type and meson-type states. Using a suitable basis of such states a matrix 
correlation can be constructed from which the static potential and its excitations are 
extracted for arbitrary separations of the static quarks. 

3.1 The matrix correlation 

The static potentials are defined as the energy levels in the charge sector of the Hilbert 
space with a static quark at space position x and a static anti-quark at space position 
x r = x + rk. States living in this charge sector are described by fields O a b(x,x r ) with 
color indices a, b and equal times xq = yo, which transform under gauge transformations 
like (HH] 

Oa b (x,x r ) = k\ al {x)O a 'b'{x,x r ) K b i h {x r ) . (22) 

The simplest choice of such fields describing string- type states is U a b(x,x r ) and for the 
meson-type states & a (x)&l(x r ). By U(x,y) we denote the product of gauge links along 
the straight line connecting y with x. The basic matrix correlation for the extraction 
of the static potentials can be expressed in terms of the following transition matrix 
elements f2§ 

C ww (r,t) = {tr[U(x,x r )U(x r ,x r + t6)U\x + tb,x r + tb)U\x,x + t6)}) , (23) 
Cwm(M) = {$Hx + tQ)tf(x,x + tO)U(x,x r )U(x r ,x r + tQ)<£>(x r + tO)), (24) 
CmmM) = (* f (a; + t6)U\x,x + t6)$(x) &(x r )U(x r ,x r + t6)$(z P +t6)) .(25) 

The static quark (anti-quark) is represented by a straight time-like Wilson line U'(x, x-\- 
£0) (U(x r ,x r + tO)). The matrix C is real, symmetric and positive. This simplest choice 
of the states does not however correspond to the physical picture that we have of the 
system. 

The string- type states should reproduce the flux tube [^7|,^8,3£,4C,41,42l of the 



gauge field binding the static quarks. Therefore, we smear the space-like links describing 



the string- type states using the APE smearing procedure of [43| with smearing strength 
set to the numerical value e = 1/4. In order to construct "physical" two-meson states, 
we determine the spectrum of the static-light mesons as described in Sect. ^. In the 
matrix correlation eq. (|5|) we use smeared Higgs fields Of L (x) = ^^"^(x) defined with 
the smearing operator of eq. (|i~l"l). The numbers rtj (i = 1, 2, N) denote the smearing 
levels. The eigenvectors v a £ 1R N (a = 0,1,2,...), obtained by solving the generalised 
eigenvalue problem eq. (0) for large t, are the wave functions describing approximately 
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(because of the finite basis of fields and the finite time t) the true eigenstates of the 
Hamiltonian. We define the fields 



N 

* a (x) = 5> Qii <4 ni) (z) (a = 0,1, 2,...), (26) 
i=i 

corresponding to the approximate one- meson eigenstates. The fields we choose to de- 
scribe two-meson states are defined as 

[* a (x)] a -[**p(x r j\ b , *,P = 0,1,2. (27) 

The values a = 0,1,2 refer to the ground, first and second excited one-meson state. 
The field basis in eq. (^) contains combinations with a ^ [3 which are not symmetric 
under interchange of the positions x and x r of the static charges. Because we expect 
the ground two-meson state to be symmetric, we project into the symmetric linear 
combinations of the fields in eq. ( p7| ) when we analyse the data of the simulations. The 
"mixed" states (for example of one meson in the ground state and one meson in the 
first excited state) can be important when looking at the asymptotic behavior (in r) of 
excited static potentials Q|. The one-meson states have a space extension due to the 
smearing of the Higgs field. For a high number of smearing iterations, there is effectively 
an "interaction" between the mesons in the two-meson states eq. (|27| ) due to the overlap 
of the smeared Higgs fields 

Summarising, we use the basis of states \i) described by the fields 

where U^ mi \x, x r ) is the product of smeared gauge links (with smearing level mi) along 
the straight line connecting x r with x and the pairs of indices (pn = 0, 1, 2; = 0, 1, 2) 
label the 9 combinations of two-meson states. We construct the following matrix corre- 
lation 

C i:j {t, r) = ([Oi{x, x r )] ab U bc {x r , x r + t6) [Oj(x + t6, x r + t0)]^ U]jx, x + t0)) . (29) 

We denote the static potentials by V a (r), a = 0, 1, 2, .... The corresponding eigenstates 
of the Hamiltonian are denoted by \a). Taking the limit of infinite time extension of 
the lattice T — ► oo, we obtain from the transfer matrix formalism the following spectral 
representation of eq. (p9j) [p^] 

Cv(t,r) = -£(j\a)(a\i)e- tv «V. (30) 

a 

For fixed separation r, we extract from C(t, r) the potentials V a (r) using the variational 
method described in Sect. P.2L 
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fields 


smearing levels m 


17 W (see ||] with e = 1/4) 


7,10,15 


(see eq. (11)) 


1,3,5,7,10,15 



Table l: Here, we list the smearing levels for the gauge and Higgs fields used in the 
simulation with parameters (3 = 2.4, k = 0.2759, A = 0.7. 



3.2 Results at /3 = 2.4 

In our first study |2^] we obtained the static potential from a simulation at (3 = 2.2, k = 
0.274, A = 0.5 on a 20 4 lattice. We observed string breaking at a distance r^/a ~ 5. 
We decided then to study the system with a better lattice resolution at j3 = 2.4. The 
results that we describe in the following are obtained on a 32 4 lattice for the parameter 
set 



P = 2.4, k = 0.2759, A = 0.7. (31) 

The field basis is constructed according to eq. (^) and eq. fl2Sj ) from smeared gauge 
(N\j = 3) and Higgs fields, whose smearing parameters are summarised in Table |l[ The 
parameters for the simulation were fixed after some trial runs. In the matrix correlation 
function Cij(t,r) = 1,2,..., 12), the time-like links are replaced by their one-link 
integrals p5fl . We collected a statistics of 800 measurements. Autocorrelations in the 
measurements are practically absent: the statistical errors, computed by a jackknife 
analysis remain constant when we group the measurements in bins of length 1,2 or 4. 



3.2.1 Renormalised static potentials 

The static potentials V a (r) (a = 0,1,2,...) are extracted from the matrix correlation 
eq. (^9|) using the variational method described in Sect. We rewrite eq. ([|) as 



aV a (r) = ln(X a (t-a,to)/X a (t,to)) + o(e- tAv ^) , (32) 

where AV a (r) = min | V a (r) — Vg(r)\ and the eigenvalues X a (t, to) are obtained by solving 

the generalised eigenvalue problem eq. (|7j) with the matrix correlation function at fixed 
r. We choose to = 0. At all distances r we can read off with confidence and very good 
statistical precision (per mille level) values for the static potential Vq{t) at t = la which 
agree fully with t = 6a. From the static potential along a lattice axis we determined 
the scale ro exactly as explained in [g]. The result is 

r /a = 5.29(6). (33) 

Comparing this number with the values of ro/a computed in quenched QCD Ej|, we 
see that our point eq. fl3l]) in the SU(2) Higgs model corresponds in resolution to (3 6 
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Figure 5: Here, the renormalised ground state and first excited state static potentials in 
units of ro are shown as functions of the separation of the static quarks. String breaking 
is clearly visible at ~ 1.9ro together with the crossing of the energy levels. 



in the SU(3) Yang-Mills theory with Wilson action. With respect to our first work [22], 
the lattice resolution is almost a factor 2 better. 



The static potentials V a (r) as they are obtained from eq. (32) are not renormalised 
quantities because they contain self-energy contributions of the static quarks which 
diverge like 1/a in the continuum. We consider instead the differences V a (r) — 2\i which 
are free of divergences |$^] and multiply them by ro to obtain dimensionless renormalised 
potentials. In Fig. ||, we represent the ground state and the first excited state static 
potentials. The ground state potential shows an approximate linear rise at intermediate 
distances: around separation 



1.9 r 



(34) 



the potential flattens: the string breaks. As expected, for large distances the potential 
approaches the asymptotic value 2fi. The first excited potential comes very close to the 
ground state potential around and rises linearly at larger distances. The scenario of 
string breaking as a level crossing phenomenon |3(J is confirmed beautifully. 
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Figure 6: In the figure on the left, we show the overlap loq obtained using the full basis 
of states eq. fl28|). In the figure on the right, we show the overlap ujff determined from 
string-type states only (described by the Wilson loops). A "naive" way of extracting 
it (triangles), using eq. (|36j), gives an erroneous large overlap at long distances. A safe 



estimate (circles) is obtained from eq. (37). 



3.2.2 Overlaps 

A certain measure for the efficiency of a basis of fields eq. (|2^) used to extract the 
ground state potential is given by the overlap. Using the approximate ground state wave 
function vq obtained from the variational method we define the projected correlation 
function 

m = E «oA-(t)«oj = £MV tv °w > ( 35 ) 

i,j a 

with normalisation f2(io = 0) = 1. The positive coefficients (w a ) 2 may be interpreted 
as the square of the overlap of the true eigenstates of the Hamiltonian \a) with the 
approximate ground state characterized by vq. The "overlap" is an abbreviation com- 
monly used to denote the ground state overlap, wq. We determine ujq straightforwardly 
from the correlation function Q(t) by noting that 

ln(u;o) 2 ~ t -^lnQ(t)--lnn(t + a). (36) 

t-*oo a a 

We extract safe values for (u>o) 2 at t = 7a, which agree fully with t = Qa and are shown 
in the left part of Fig. |(| Our basis of fields is big (and good) enough such that (o>o) 2 
exceeds about 60% for all distances. 

It is interesting to consider also the overlap for the (smeared) Wilson loops alone, 
i.e. we restrict the matrix correlation function to the 3x3 sub-block associated with 
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Figure 7: In the figure on the left, the overlaps of the string-type (circles) and meson- 
type (triangles) states with the ground state of the Hamiltonian are shown as functions 
of the separation r of the static quarks. In the figure on the right, the same overlaps 
but with the first excited eigenstate of the Hamiltonian are shown. 



string-type states. Let us denote the corresponding projected correlation function by 
OwW and the overlap by off. The computation of uiff is more difficult because it turns 
out to be very small at large r. In the right part of Fig. ^, we present the results for 
two estimates of (uSq) 2 . The triangles correspond to the estimates (w^ V,naive ) 2 obtained 
directly from eq. (|36|), with O(t) replaced by Qw(t). The circles correspond to the more 
reliable estimate using the information from the full matrix correlation: the expression 

K w ) 2 ~ M 2 %^ (37) 

converges reasonably fast and (uj^) 2 can be estimated from the r.h.s. for large t (t/a = 
7 — 9 in practice). Using eq. (|37|), we see that (smeared) Wilson loops alone have an 
overlap which drops at intermediate distances and they are clearly inadequate to extract 
the ground state at large r. On the contrary, (o;J V ' naivc ) 2 is above 50% at large distances: 
what is estimated here, is actually the coefficient (u>^) 2 , i.e. the square of the overlap 
of the (smeared) Wilson loops with the first excited state (this statement is supported 
by direct calculation, see Sect. P-2.3Q . The fact that uj^ is so large might explain 



the problems encountered in QCD for observing string breaking from the analysis of a 
correlation matrix with Wilson loops only. 

3.2.3 Mixing 

Finally, we want to show that string breaking is a mixing phenomenon involving string- 
type and meson-type states. This leads to the crossing of the energy levels seen above. 
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We consider the diagonal sub-blocks of the matrix correlation function eq. ( |29|) 
corresponding to string- type states (fields i = 1,2,3 in eq. (p8[)) and to meson- type 
states (fields i = 4, 5,. ..,12 in eq. (p8|)) separately. We determine approximate ground 
state wave functions vff in the subspace of the string-type states and Vq 1 in the subspace 
of the meson-type states. With the help of these wave functions we construct a 2 x 2 
projected matrix correlation function 

= E «0,iCtf(*Kj = E^H(a|^)e-^ (r) (M = W,M) , (38) 

where 

Mk> = 5>o,il*>- (39) 

i 

An inspection of 

n WM (t = o) = (VmI^w) (40) 

shows that string-type and meson-type states are orthogonal only for large values of r 
. The coefficients 

w fc (a) = (a\^ k ) (fc = W,M), (41) 

in the expansion eq. (|38|), express the overlap of the string-type (k = W) and meson- 
type (k = M) states with the true eigenstates of the Hamiltonian. We can choose 
our phase conventions for the states such that the coefficients u)k{a) are real and in 
addition i%(0) > and i%(l) > 0. We truncate the sum in eq. (p8|) after a = 1 and 
consider the diagonal matrix elements Slkk{t) for two fixed times t = t\ and t = t2- 
inserting the known values for Vo(r) and Vi(r), we can solve for w|(0) and cc?(l). The 
sign of the coefficients <^m(0) and wm(1) is fixed by the off-diagonal matrix elements 
^wm(^i) and Owm(*2): we find that for all r, wm(0) > and wm(1) < (in our sign 
convention). The overlaps u>-w(0) of the string- type states (circles) and wm(0) of the 
meson-type states (triangles) with the ground state of the Hamiltonian are shown on 
the left of Fig. and the corresponding overlaps and wm(1) with the first excited 

eigenstate of the Hamiltonian on the right of Fig. ^. String-type states have a large 
overlap at short distances with the ground state and at large distances with the first 
excited state. Meson-type states have a large overlap at short distances with the first 
excited state and at large distances with the ground state. 

In addition, we observe that the overlap of the meson- type states with the ground 
state is also large at very short distances. The explanation for this fact is that string- 
type and meson-type states have an overlap with each other at short distances. In the 
string breaking region around r/a = 9 — 10, the overlaps of the string-type and meson- 
type states have similar magnitude, both when the ground state or the first excited state 
is considered. This fact is reflected in the crossing of the energy levels Fig. ||. Here, we 
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Figure 8: Here, we show the scaling of the renormalised ground state and first excited 
state static potentials. The parameter points lie on a line of constant physics. 

would like to point out that the overlaps represented in Fig. ^ are not quantities which 
have a strict continuum limit. They depend on the (5- value and the other parameters 
(e.g. of the smearing) that we consider. However, as long as one chooses a good basis 
(say with luq > 0.5) which can be separated into "string like" and "meson like", the 
qualitative behavior in Fig. |7| is expected to persist also at smaller lattice spacings. 

4 Scaling of the static potentials 

In order to compare the renormalised static potentials at different values of the lattice 
spacing and estimate the size of scaling violations, a way of determining lines of constant 
physics (LCP) in the confinement "phase" of the SU(2) Higgs model is needed. This 
question has been addressed in [31], where a non-perturbative determination of the 
LCPs is described. 

The bare parameters k and A are renormalised along a LCP by keeping two physical 
quantities F\ and F2 constant. A good choice is to take 



F x = r [2/x - Vb(ro)] 
and F2 to be the generalised Binder cumulant C3 defined in [^IJ. f\ 



(42) 



7 Physically Fx is a (non-perturbative) measure of the Higgs mass: the variation of Fx with the 
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Figure 9: Here, we show the A-independence of the static potentials. The parameter k 
is determined from eq. ([45|). 



Taking the parameter set (^Tj) which we used in the simulation at /? = 2.4 we obtain 
the conditions 



F* 



1.26 and F 2 



(43) 

As a result of the 



defining a LCP (the numerical value F£ can be found in plf ). 
non-perturbative matching, the parameter sets ([H]) and 

/? = 2.2, k* = k(X*), A* = 0.96(10) (44) 

lie on a LCP. The value of k at (3 = 2.2 is determined using the polynomial fit |3l| 

k(A) = 0.3131 + 0.0564 (A - 1) - 0.0286 (A - l) 2 

+ 0.0198 (A - l) 3 -0.0246 (A - l) 4 (45) 

which is obtained by requiring F% = F* = 1.26 and correlates the uncertainties in k and 
A of the LCP. In Fig. || we compare the results for the renormalised ground state and 

bare parameters is dominantly caused by the variation of the meson mass {i. This bound state mass 
is of course expected to depend strongly on the mass of its constituents. On the other hand, the 
interpretation of F2 is less obvious. It was chosen to have a second renormalised quantity which is both 
sensitive to the bare coupling A and can be computed in the MC simulations fell. 
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first excited state static potentials that we obtain at (5 = 2.4 and at (3 = 2.2 along the 
LCP. The results for the potentials are compatible with scaling within minute errors 
under variation of the lattice spacing by almost a factor 2. 

Another interesting issue is the A-dependence of physical observables. Exploratory 
results |Q indicated that the physics of the SU(2) Higgs model in the confinement 
"phase" is weakly dependent on A. This is certainly true near the continuum limit 
because it is well accepted that the scalar part of the SU(2) Higgs model is a trivial 
theory. Nevertheless, at finite value of the lattice spacing the model can be considered 
as an effective field theory with three independent renormalised couplings. In Fig. [9], we 
compare the static potentials for two different values A = 0.5 and A = 0.96 at [3 = 2.2. 
The parameter k is determined from eq. (f45|). There is no significant difference between 
the different A values. This also means that the uncertainty in A* in ( (44] ) is irrelevant 
for our scaling test Fig. ^. 

5 Conclusions 

We presented the results for the static potentials (ground state and first excited state) 
in the confinement "phase" of the SU(2) Higgs model. The string breaking or flattening 
of the ground state potential is clearly visible around separation ~ 1.9ro. The 
comparison with the first excited potential shows a nice crossing of the energy levels. 
The interpretation of string breaking as level crossing phenomenon between string-type 
and meson-type states is substantiated by the investigation of properly defined overlaps. 

We also addressed the question of scaling violations in the measurements of the 
static potentials. They are shown to be tiny already at r^ja > 2.5. Moreover, the 
dependence of the static potentials on the Higgs quartic coupling A is very weak once 
the parameter k is determined by keeping the physical quantity iq constant. These 
results are a strong indication for a continuum-like behavior of the static potentials 
already at the relatively large values of the lattice spacing that we used. 

The method for the determination of the static potential that we presented can be 
applied in QCD. Some steps in this direction have already been made [p8| . The main 
problem for QCD is the statistical accuracy. 

Acknowledgement. We thank the Konrad-Zuse-Zentrum fur Informationstechnik 
Berlin (ZIB) for granting CPU-resources to this project. 
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